The observed velocity distribution of young pulsars – II. Analysis of complete PSRπ

Understanding the natal kicks, or birth velocities, of neutron stars is essential for understanding the evolution of massive binaries and double neutron star formation. We use maximum likelihood methods as published in Verbunt et al. to analyse a new large data set of parallaxes and proper motions measured by Deller et al. This sample is roughly three times larger than number of measurements available before. For both the complete sample and its younger part (spin-down ages τ < 3 Myr), we find that a bimodal Maxwellian distribution describes the measured parallaxes and proper motions better than a single Maxwellian with probability of 99.3 and 95.0 per cent, respectively. The bimodal Maxwellian distribution has three parameters: fraction of low-velocity pulsars and distribution parameters σ1 and σ2 for low- and high-velocity modes. For a complete sample, these parameters are as follows: 42+17−15 per cent, σ1=128+22−18 km s−1, and σ2 = 298 ± 28 km s−1. For younger pulsars, which are assumed to represent the natal kick, these parameters are as follows: 20+11−10 per cent, σ1=56+25−15 km s−1, and σ2 = 336 ± 45 km s−1. In the young population, 5 ± 3 per cent of pulsars have velocities less than 60 km s−1. We perform multiple Monte Carlo tests for the method taking into account realistic observational selection. We find that the method reliably estimates all parameters of the natal kick distribution. Results of the velocity analysis are weakly sensitive to the exact values of scale lengths of the Galactic pulsar distribution

Braking indices of young radio pulsars: theoretical perspective

Recently, Parthsarathy et al. analysed long-term timing observations of 85 young radio pulsars. They found that 15 objects have absolute values of braking indices ranging ∼10 − 3000, far from the classical value n = 3. They also noted a mild correlation between measured value of n and characteristic age of a radio pulsar. In this article we systematically analyse possible physical origin of large braking indices. We find that a small fraction of these measurements could be caused by gravitational acceleration from an unseen ultra-wide companion of a pulsar or by precession. Remaining braking indices cannot be explained neither by pulsar obliquity angle evolution, nor by complex high-order multipole structure of the poloidal magnetic field. The most plausible explanation is a decay of the poloidal dipole magnetic field which operates on a time scale ∼104 − 105 years in some young objects, but has significantly longer time scale in other radio pulsars. This decay can explain both amplitude of measured n and some correlation between n and characteristic age. The decay can be caused by either enhanced crystal impurities in the crust of some isolated radio pulsars, or more likely, by enhanced resistivity related to electron scattering off phonons due to slow cooling of low-mass neutron stars. If this effect is indeed the main cause of the rapid magnetic field decay manifesting as large braking indices, we predict that pulsars with large braking indices are hotter in comparison to those with n ≈ 3

Inferred time-scales for common envelope ejection using wide astrometric companions

Evolution of close binaries often proceeds through the common envelope stage. The physics of the envelope ejection (CEE) is not yet understood, and several mechanisms were suggested to be involved. These could give rise to different time-scales for the CEE mass-loss. In order to probe the CEE-time-scales we study wide companions to post-CE binaries. Faster mass-loss time-scales give rise to higher disruption rates of wide binaries and result in larger average separations. We make use of data from Gaia DR2 to search for ultrawide companions (projected separations 103–2 × 105 au and M2 > 0.4 M⊙) to several types of post-CEE systems, including sdBs, white dwarf post-common binaries, and cataclysmic variables. We find a (wide-orbit) multiplicity fraction of 1.4 ± 0.2 per cent for sdBs to be compared with a multiplicity fraction of 5.0 ± 0.2 per cent for late-B/A/F stars which are possible sdB progenitors. The distribution of projected separations of ultrawide pairs to main sequence stars and sdBs differs significantly and is compatible with prompt mass-loss (upper limit on common envelope ejection time-scales of 102 yr). The smaller statistics of ultrawide companions to cataclysmic variables and post-CEE binaries provide weaker constraints. Nevertheless, the survival rate of ultrawide pairs to the cataclysmic variables suggest much longer, ∼104 yr time-scales for the CEE in these systems, possibly suggesting non-dynamical CEE in this regime

Powering Central Compact Objects with a Tangled Crustal Magnetic Field

Central Compact Objects (CCOs) are X-ray sources with luminosity ranging between 1032-1034 erg s−1, located at the centres of supernova remnants. Some of them have been confirmed to be neutron stars. Timing observations have allowed the estimation of their dipole magnetic field, placing them in the range ∼1010-1011 G. The decay of their weak dipole fields, mediated by the Hall effect and Ohmic dissipation, cannot provide sufficient thermal energy to power their X-ray luminosity, as opposed to magnetars whose X-ray luminosities are comparable. Motivated by the question of producing high X-ray power through magnetic field decay while maintaining a weak dipole field, we explore the evolution of a crustal magnetic field that does not consist of an ordered axisymmetric structure, but rather comprises a tangled configuration. This can be the outcome of a non-self-excited dynamo, buried inside the crust by fallback material following the supernova explosion. We find that such initial conditions lead to the emergence of the magnetic field from the surface of the star and the formation of a dipolar magnetic field component. An internal tangled magnetic field of the order of 1014 G can provide sufficient Ohmic heating to the crust and power CCOs, while the dipole field it forms is approximately 1010 G, as observed in CCOs

Hyper-runaway and hypervelocity white dwarf candidates in Gaia Data Release 3: possible remnants from Ia/Iax supernova explosions or dynamical encounters

Type Ia and other peculiar supernovae (SNe) are thought to originate from the thermonuclear explosions of white dwarfs (WDs). Some of the proposed channels involve the ejection of a partly exploded WD (e.g. Iax SN remnant) or the companion of an exploding WD at extremely high velocities (>400 km s−1). Characterisation of such hyper-runaway/hypervelocity (HVS) WDs might therefore shed light on the physics and origins of SNe. Here we analyse the Gaia DR3 data to search for HVS WDs candidates, and peculiar sub-main-sequence (sub-MS) objects. We retrieve the previously identified HVSs, and find 46 new HVS candidates. Among these we identify two new unbound WDs and two new unbound sub-MS candidates. The remaining stars are hyper-runaway WDs and hyper-runaway sub-MS stars. The numbers and properties of the HVS WD and sub-MS candidates suggest that extreme velocity ejections (>1000 km s−1) can accompany at most a small fraction of type Ia SNe, disfavouring a significant contribution of the D6-scenario to the origin of Ia SNe. The rate of HVS ejections following the hybrid WD reverse-detonation channel could be consistent with the identified HVSs. The numbers of lower-velocity HVS WDs could be consistent with type Iax SNe origin and/or contribution from dynamical encounters. We also searched for HVS WDs related to known SN remnants, but identified only one such candidate

Strong toroidal magnetic fields required by quiescent X-ray emission of magnetars

Magnetars are neutron stars (NSs) with extreme magnetic fields1 of strength 5 × 1013−1015 G. These fields are generated by dynamo action during the proto-NS phase, and are expected to have both poloidal and toroidal components2,3,4,5,6, although the energy of the toroidal component could be ten times larger7. Only the poloidal dipolar field can be measured directly, via NS spin-down8. The magnetic field provides heating and governs how this heat flows through the crust9. Magnetar thermal X-ray emission in quiescence is modulated with the rotational period of the NS, with a typical pulsed fraction 10–58%, implying that the surface temperature is substantially non-uniform despite the high thermal conductivity of the star’s crust. Poloidal dipolar fields cannot explain this large pulsed fraction10,11. Previous two-dimensional simulations12,13 have shown that a strong, large-scale toroidal magnetic field pushes a hot region into one hemisphere and increases the pulsed fraction. Here, we report three-dimensional magneto-thermal simulations of magnetars with strong, large-scale toroidal magnetic fields. These models, combined with ray propagation in curved spacetime, accurately describe the observed light curves of 10 out of 19 magnetars in quiescence and allow us to further constrain their rotational orientation. We find that the presence of a strong toroidal magnetic field is enough to explain the strong modulation of thermal X-ray emission in quiescence

Combined analysis of neutron star natal kicks using proper motions and parallax measurements for radio pulsars and Be X-ray binaries

Supernova explosion and the associated neutron star (NS) natal kicks are important events on a pathway of a binary to become a gravitational wave source, an X-ray binary, or a millisecond radio pulsar. Weak natal kicks often lead to binary survival, while strong kicks frequently disrupt the binary. In this article, we aim to further constrain NS natal kicks in binaries. We explore binary population synthesis models by varying prescription for natal kick, remnant mass, and mass accretion efficiency. We introduce a robust statistical technique to analyse combined observations of different nature. Using this technique, we further test different models using parallax and proper motion measurements for young isolated radio pulsars and similar measurements for Galactic Be X-ray binaries (BeXs). Our best model for natal kicks is consistent with both measurements and contains a fraction of w = 0.2 ± 0.1 weak natal kicks with σ1=45+25−15 km s−1, the remaining natal kicks are drawn from the high-velocity component, same as in previous works: σ2 = 336 km s−1. We found that currently used models for natal kicks of NSs produced by electron capture supernova (ecSN; combination of Maxwellian σ = 265 km s−1 and σ = 30 km s−1 for electron capture) are inconsistent or marginally consistent with parallaxes and proper motions measured for isolated radio pulsars. We suggest a new model for natal kicks of ecSN, which satisfy both observations of isolated radio pulsars and BeXs

3D Magnetothermal Simulations of Tangled Crustal Magnetic Field in Central Compact Objects

Central compact objects (CCOs) are young neutron stars emitting thermal X-rays with bolometric luminosities LX in the range of 1032–1034 erg s−1. Gourgouliatos, Hollerbach, and Igoshev recently suggested that peculiar emission properties of CCOs can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto–neutron star stage. In this case the magnetic field consists of multiple small-scale components with negligible contribution of global dipolar field. We study numerically three-dimensional magnetothermal evolution of tangled crustal magnetic fields in neutron stars. We find that all configurations produce complicated surface thermal patterns that consist of multiple small hot regions located at significant separations from each other. The configurations with initial magnetic energy of (2.5–10) × 1047 erg have temperatures of hot regions that reach ≈ 0.2 keV, to be compared with the bulk temperature of ≈ 0.1 keV in our simulations with no cooling. A factor of two in temperature is also seen in observations of CCOs. The hot spots produce periodic modulations in light curve with typical amplitudes of ≤9%–11%. Therefore, the tangled magnetic field configuration can explain thermal emission properties of some CCOs

Discovery of X-Rays from the Old and Faint Pulsar J1154-6250

We report on the first X-ray observation of the 0.28 s isolated radio pulsar PSR J1154-6250 obtained with the XMM-Newton observatory in 2018 February. A point-like source is firmly detected at a position consistent with that of PSR J1154-6250. The two closest stars are outside the 3 sigma confidence limits of the source position and thus unlikely to be responsible for the observed X-ray emission. The energy spectrum of the source can be fitted equally well either with an absorbed power law with a steep photon index Gamma approximate to 3.3 or with an absorbed blackbody with temperature kT = 0.21 +/- 0.04 keV and emitting radius R-BB approximate to 80 m (assuming a distance of 1.36 kpc). The X-ray luminosity of 4.4 x 10(30) erg s(-1) derived with the power-law fit corresponds to an efficiency of eta(X) = L-X(unabs) /(E) over dot= 4.5 x 10(-3), similar to those of other old pulsars. The X-ray properties of PSR J1154-6250 are consistent with an old age and suggest that the spatial coincidence of this pulsar with the OB association Cm OB1 is due to a chance alignment